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Abstract 

We compute a formula including OPE power corrections to describe the running of a 
QCD coupling non-perturbatively defined through the ghost and gluon dressing functions. 
This turns out to be rather accurate. We propose the "p^aieau" -procedure to compute A^g 
from the lattice computation of the running coupling constant. We show a good agreement 
between the different methods which have been used to estimate • We argue that 

MS to 

Aj^g or the strong coupling constant computed with different lattice spacings may be used 
to estimate the lattice spacing ratio. 

LPT-Orsay 08-91 
UHU-FT/08-10 

1 Introduction 

Much work has been devoted in the last years to the study of the QCD running coupling con- 
stant determined from lattice simulations, as well in its perturbative regime [1-9] as in the deep 
infrared domain [10]. The two main approaches to obtain the running coupling in terms of the 
renormalization momentum were either an application of the Schrodinger functional method with 
special boundary conditions or the confrontation of the behaviour with respect to the renormaliza- 
tion scale of 2-gluon and 3-gluon Green functions with the corresponding perturbative predictions. 
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The latter Green's functions approach also revealed a dimension-two non-zero gluon condensate in 
the landau gauge. Much work has been also done to investigate its phenomenological implications 
in the gauge- invariant world [11]. In a very recent work [9], the Green's function approach to 
estimate A^g has been pursued by exploiting a non-perturbative definition of the coupling derived 
from the ghost-gluon vertex and computed over a large momentum window in the perturbative 
regime. Much of this work was based on the analysis of quenched lattice simulations and led to the 
determination of A^jg in pure Yang-Mills (Nf = 0). Works on unquenched lattice configurations 
(Nf = 2) started some time ago [12] and have been more actively pursued recently. 

Many unquenched configurations are presently available and we are planing to apply what we 
have learned on pure Yang-Mills to gauge configurations with twisted Nf = 2 [13]) and Nf = 
2 + 1 + 1 dynamical quarks. Thus, a very realistic estimate of A^g-, directly comparable with 
experimental determinations, will become an immediate possibility. With the latter remarks in 
mind, we pay attention in this paper to study the above-mentioned non-perturbative coupling 
derived from the ghost-gluon vertex for being applied to the analysis of lattice data. We show in 
section 2 that, when the incoming ghost-momentum vanishes -and only in this case- this ghost- 
gluon vertex can be directly related to the bare gluon and ghost propagators; we then obtain 
a formula to describe its running including non-perturbative power corrections. We propose to 
confront this formula with lattice estimates of the coupling and argue that this constitutes an 
optimal method for the identification of A^jg and of the gluon condensate. In particular, it 
benefits of two main advantages: to have only two-points function to deal with (much simpler to 
be managed and more precise than three-points ones) and that the precision could be improved by 
extending the analysis of lattice data over a very large momenta window. In section 3, we apply 
this procedure to previouly published lattice data for quenched simulations with a two-sided goal: 
(i) to check the method and (ii) to confirm the consistency of the picture we have acquired for the 
UV behaviour of Green functions in pure Yang-Mills. We finally conclude in section 4. 

2 The ghost-gluon coupling 

There is a large number of possibilities to define the QCD renormalized coupling constant, de- 
pending on the observable used to measure it and on the renormalization scheme. Actually, any 
observable which behaves, from the perturbative point of view, as g provides a suitable defini- 
tion for it. Among such quantities stand the 3-gluon and the ghost-gluon vertices, which have 
been widely used by the lattice community to get a direct knowledge of a s from simulations. Of 
course an important criterion to choose among those definitions will be how easy it is to connect 
it to other commonly used definitions, specially the MS one, and to extract from it fundamental 
parameters like Aq C d- 

A convenient class of renormalization schemes to work with on the lattice is made of the so- 
called "MOM" schemes which are defined through the requirement that a given scalar coefficient 
function of the Green's function under consideration take its tree-level value in a specific kinemat- 
ical situation given up to an overall "renormalization scale" . To make the point clearer we recall 
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2 schemes which we have used in previous works on a s : 

• The symmetric 3-gluon scheme in which one uses the 3-gluon vertex T tlvp (pi,p2,P3) with 

v\ = vl = vl = 

• The asymmetric 3-gluon scheme (MOM) in which the 3-gluon vertex ^ P up(PiiP2,P3) is used 
with p\ = p\ = fi 2 , p\ = 

In the present note we shall apply a specific MOM-type renormalization scheme defined by 
fixing the (ghost and gluon) propagators and the ghost-gluon vertex at the renormalization point. 
Let us start by writing the ghost and gluon propagators in Landau gauge as follows, 

(GW)>>) = , 

(F<TV.A) = -4*^^; (l) 

A being some regular isat ion parameter (a _1 (/3) if, for instance, we specialise to lattice regularisa- 
tion). The renormalized dressing functions, G R and F R are defined through : 



G R {p 2 ^ 2 ) = \im Z^ 2 , A) G(p 2 , A) 

A-^oo 

F R (p 2 ^ 2 ) = lim Z^(^,A) F(p 2 ,A) , (2) 

A— >oo 

with renormalization condition 

G R {n 2 ,n 2 )=F R {n 2 ,n 2 ) = l. (3) 

Now, we will consider the ghost-gluon vertex which could be non-perturbatively obtained through 
a three-point Green function, defined by two ghost and one gluon fields, with amputated legs after 
dividing by two ghost and one gluon propagators. This vertex can be written quite generally as: 



i q_k 

V fc ( k . q - k)= _^5^_ = igof abc ( fc) + ( _ k ) vIl2 ( q k) ) ( 4 ) 

k q 

where q is the outgoing ghost momentum and k the incoming one, and renormalized according to: 



Z 1 T. 



(5) 



The vertex T y involves two independent scalar functions. In the MOM renormalization procedure 
Zi is fully determined by demanding that one specific combination of those two form factors 
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(chosen at one's will) be equal to its tree-level value for a specific kinematical configuration. We 
choose to apply MOM prescription for the scalar function Hi + H 2 that multiplies q v in eq. (j4|) 
and the renormalization condition reads 2 ! 



(H*(q,k) + HK(q,k))\ = lim Z^\K) (Hi(q, k; A) + H 2 (q, k; A))| 2 2 = 1, (6) 

where we prescribe a kinematics for the substraction point such that the outgoing ghost momentum 
is evaluated at the renormalization scale, while the incoming one, k, depends on the choice of 
several possible configurations; for instance: k 2 = (q — k) 2 = fi 2 (symmetric configuration) or 
k — 0, (q — k) 2 = fi 2 (asymmetric-ghost configuration). 

On the other hand, the fields involved in the non-perturbative definition of the vertex Y u in 
eq. (j3J) can be directly renormalized by their renormalization constants, Z 3 and Z 3 , and the same 
MOM prescription applied to the scalar combination H\ + H 2 also implies: 

g R (fi 2 ) = lim Z 3 ^ 2 ,A)Z 1 3 /2 ( f i 2 ,A)g (A 2 ) (H 1 (q,k;A) + H 2 (q,k;A) 

A— >oo \ 

= i im go( A*) gi^SM . ' ' (7) 

Zi((i 2 ,A 2 ) 

We combine both eq. (El) and the first-line equation of (JTj) to replace Hi + H 2 and obtain the 
second line that shows the well-known relationship Z g = {Z\ Zz) -1 Zi, where = Z^go- 

We turn now to the specific MOM-type renormalization scheme defined by a zero incoming 
ghost momentum. Since those kinematics are the ones (and the only ones) in which Taylor's 
well known non- renormalization theorem (cf. ref [14]) is valid we shall refer to this scheme as to 
the T-scheme and the corresponding quantities will bear a T subscript. Then, in eq 01]), we set k 
to and get 

K bc (-q, 0; q) = tg f abc (#i(g, 0) + H 2 (q, 0)) q v . (8) 

Now, Taylor's theorem states that Hi(q, 0; A) + H 2 (q, 0; A) is equal to 1 in full QCD for any value 
of q. Therefore, the renormalization condition eq. (EJ) implies Zi(fi 2 ) = 1 and then 

artf) ee 9 ^fl = nm g -Ml G (n 2 , A 2 )F V, A 2 ) ; (9) 
Air A->oo Air 



where we also apply the renormalization condition for the propagators, eqs. ( 121I31) . to replace the 
renormalization constants, Z 3 and Z3, by the bare dressing functions. The remarkable feature of 
eq. (Q is that it involves only F and G so that no measure of the ghost-gluon vertex is needed for 
the determination of the coupling constant. 



2 In the case of zero-momentum gluon, an appropriate choice would be Zi([i 2 )Hi(q,q)\ q 2 = ^2 = 1. This would 
make the renormalized vertex equal to its tree-level value at the renormalization scale. 
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Equation ([9]) has extensively been advocated and studied on the lattice (see for instance refer- 
ence [15]) and used for a determination of Aqcd in reference [9]. However it must be stressed that 
the T-scheme is the only one in which Z\ = 1. In any other scheme Z\ will be finite (since going 
from one scheme to any other one only involves an additional finite renormalization) but will keep 
a non trivial dependence on the scale, in particular for the symmetric scheme of reference [16] 
that has been computed at one loop in ref. [17]. In such cases one must in principle apply the 
general definition (J7|) of the coupling constant; nevertheless the form @ is used quite often in this 
case (for a kinematical configuration other than T-scheme's) also as an approximation, specially 
in relation with the study of Dyson-Schwinger equations. 

We conclude this section by recalling that, in any scheme, the standard renormalization flow 
dictating the evolution with respect to the scale, 

gW) = gW 2 ) (§^) ^(/A/O^V 2 ), (io) 

will be straightfowrdly obtained from the second line of eq. (J7]) and the propagators renormalization 
conditions in eqs. ( I2l3l) . where 

Zx(/i 2 ) = lim Zx(/i 2 ,A 2 ) (11) 

A^oo 

because of the Taylor's non-renormalization theorem. Of course, eq. f llOp reduces to 

gW) = gW 2 )F 2 R (» 2 ,»' 2 )G R (» 2 ,»' 2 ). (12) 

in the T-scheme. 



2.1 Pure perturbation theory 

In ref. [18], the three-loop perturbative substraction of all the three- vertices appearing in the 
QCD Lagrangian for kinematical configurations with one vanishing momentum has been done (in 
particular, the one involved in the definition of the coupling by eq. 0). Different definitions of 
the coupling constant can be related in perturbation theory through relations like : 

a T (^) = a(^) (l + ^cJ^p-X) ; (13) 



8=1 



on the other hand, since eq. ([9]) completely defines the running of the coupling, after properly 
deriving both its l.h.s. and r.h.s., one obtains 

1 daxiu 2 ) 1 /„ , d , „. 9 d 



Q!t([J> 2 ) da /^Ms( a ) V A-^oorfln/i 2 ' A^oorfln/i 2 

27(a) + 7(a) 



2 lim -p^lnF^A) + hm -j^--^ In G^ 2 , A) 

(14) 
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where 



da \ ^ I a 



= 5w? = - 4ir 5 ft UJ (15) 



is the standard /3-function for the running coupling renormalized according to the usual MS 
prescription, while 

^lnZ 3;M OM(/i 2 , A) r (ilnF(/i 2 , A) 
7W ~ A™ dln//2 ~ A™ dln/i^ _ 

c?lnZ 3M OM(^ 2 , A) rflnG(^ 2 ,A) ^ / a V +1 

71a) = hm = hm = — > 7; — 16 

,v 7 A-oo rfln/i 2 A-,00 dln/1 2 V 47r / 

are the anomalous dimensions for gluon and ghost propagators, both renormalized along MOM 
prescriptions (i.e., G^(/i 2 ,/i 2 ) = Fft(/i 2 ,/i 2 ) = 1), but expanded in terms of the MS coupling a. 
The (5 i coefficients in eq. (fl5j) have been computed up to four loops in ref. [19], (3 Q and (3 X being 
scheme- independent. Then, eqs. f)13|15|[TB|) can be applied to eq. (THl) and one is led to deal with 
a coupled system of n algebraic equations to compute the coefficients q and determine ax at 
n loops. To summarize, the running of coupling constant ax, although formally defined from a 
three-point Green function, can be derived from the knowledge of the standard MS /3-function and 
only two-points functions for ghost and gluon. These two anomalous dimensions were computed 
in the MS scheme at four loops in ref. [20] and were converted into the MOM scheme in ref. [21] 
for Nf = by applying 

r dln(Z rMS (fi 2 ,A)) din {T Rm (a)) 

lv,MOM{a) = hm — + — 

a^oo a in n z a In /i z 

= 7 « ( " ) + rfl^^ lnr » ( " ) ' (17) 

where T stands generically for the two bare two-point dressing functions F and G, for the 
renormalized ones 1 and Z r for the appropriate ^normalization constant. Eq. (E) provides also 
the coefficients 7, and 7, for any Nf (see appendix [X]) . Thus, one can solve the above mentioned 
coupled system of algebraic equations and obtain the coefficients q in eq. (fT3|) . the first of those 
equations (the one stemming from matching the 1 / a-terms in the 2 sides) resulting in the following 
constraint □: 

27o + 7o = ^o, (18) 



3 The gluon and ghost renormalized propagators in the MS scheme were also provided by ref. [18] 
4 Eq. psp is a well-known relation verified by scheme-independent coefficients of the ghost and gluon anomalous 
dimensions and of the /^-function. 
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which, in this context, results from eq. ([9]). The three first coefficients Cj in Landau gauge, for 
instance, will be: 



507 - A0N f 
Cl " 36 



76063 351 , , /1913 4 , A Ar 100 „ r9 

c 2 = C 3 - + -C 3 N f + N f 

2 144 8 SV ; V 27 3 / 91 

42074947 60675 A , o . 70245 . _ Z769387 8362 A . Q . 2320 _, 

C3 = "^728 16- C(3) --64- C(5) - ^2--TT C(3) -^ C(5)) N ' 



+ ^7- + irC(3) (19) 



/ 199903 28. /nN \ Ar2 1000 , 

1 C(3) iV? iV? 

V 972 9 SV 7 J 729 f 

These three coefficients obviously define unambigously the running of ay given in eq. ([9]) up to 
four-loops. In other words, one obtains for the /^-function of ar, 



i=0 

the following coefficients up to four-loops 
A, = P = ll-^N f 

A = ^ = 102-^ 

Aa = 73 2 -7Vi + A)(c2-Ci) 

= 3040.48 - 625.387 N f + 19.3833 iVj 

& = 7i 3 -2^ 2 Ci+7iiC? + 7io(2 c 3 -6 c 2Cl + 4 C 3) 

= 100541 - 24423.3 N f + 1625.4 iVj - 27.493 Nf , (21) 

These coefficients $ are the same as the ones obtained in ref. [18] thanks to a direct application 
of the MOM prescription to the ghost-gluon coupling with vanishing incoming-ghost momentum, 
as it should be. As for the Aqcd parameters in the two schemes, they are related through 

ci 507 - A0N f 

%s = = e ~792 -48JV / (22) 

At 
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Eq. ( !20l) can be integrated and perturbatively inverted to obtain the following standard four-loop 
formula for the running coupling: 

2 4n ( ft log(i) Bf 1 / / lV.AAi 5\ \ 

+ us? (4 + Kl) 3 (- 2 ,og3w + 5 iog2(() + ( 4 - 6 f ) ,ogW - i ) ) (23) 

with t = In -^r . 

Af, 

As a last remark, applying the approximation Z\ — 1 for symmetric (ghost-gluon vertex 
renormalized at a symmetric momenta configuration) or soft-gluon (vertex renormalized at a 
vanishing-gluon momenta configuration) schemes implies that the same lattice data for the cou- 
pling, obtained through eq. @, would be confronted to different perturbative formulae analogous 
to eq. ( 123*1) with /3-function coefficients and Aqcd parameters apropriate for each scheme. Thus, 
the systematic deviation induced by applying this approximation to the determination of A^jg from 
the confrontation of perturbation theory and lattice data, provided that (3q and f3\ are scheme- 
independent, mainly results from the ratio of Aq CD to A^jg in eq. ([""""") . For instance in pure 
Yang-Mills, if one takes Nf = in eq. (1221) . it gives a ratio of 0.527 in T-scheme, while the same 
ratio for instance in symmetric and soft-gluon schemes is 0.463 (14 % of error) and 0.429 (23 % 
of error), respectively. 

2.2 OPE power corrections 

One of the goals of the present paper consists in obtaining a formula for the QCD running coupling 
that could be implemented in conjunction with lattice estimates to determine a "plateau" for 
Aqcd in terms of the momentum, as will be explained in the next section. In order to extend this 
"plateau" to energies as low as possible (of the order of 3 GeV) and to take full advantage of the 
lattice data in order to reduce the systematic uncertainties, it is mandatory to take into account 
the gauge- dependent dimension-two OPE power corrections (cf. [7,8,10,23]) to ay. 
The leading power contribution to the ghost propagator, 

(F^) ab (q 2 ) = J <$xe iq x { T (c a (x)?(0)) ) (24) 

can be computed using the operator product expansion [24] (OPE), as is done in ref. [22], 

T(c a (x)?(0))=^( Ci r b (x)O t (0); (25) 

t 

here O t is a local operator, regular when x — > 0, and the Wilson coefficient c t contains the 
short-distance singularity. Eq. (]25p involves a full hierarchy of terms, ordered according to their 
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mass-dimension, among which only 1 and : A^A b u : contribute to eq. ( l24j) in Landau gauge up 
to the order l/q 4 . Then, using eq. (1251) into eq. (El]), we obtain: 



(^ (2) ) a V) = (co) a V) + (c 2 )7 T (<? 2 )<: A'(0)4(0):> 



(iffiW) + ^ TTl^TT + • • • (26) 



where 



i / dW** (A*',(0) T (c a c*>) ^',(0)) comiected 



2 G( 2 ) ss ',G( 2 ) tt ', 

CTCr' TT 



= 2x g„_ , (27) 

and the SVZ factorisation [25] is invoked to compute the Wilson coefficients. Thus, one should 
compute the "sunset' diagram of the last line of eq. (j2*7|) . that binds the ghost propagator to the 
gluon condensate (where the blue bubble means contracting the color and lorentz indices of the 
incoming legs with l/2S s tS aT ) to obtain the leading non-perturbative contribution (of course, the 
first Wilson coefficient gives trivially the perturbative propagator). 
Finally, 



(^rW) = (i^ ert rv^ 2 ) (i+^ f^g ) + (^" 4 ) (28) 

where the A 2 -condensate is renormalised at the subtraction point q 2 = fi 2 , according to the MOM 
scheme definition, by imposing the tree-level value to the Wilson coefficient at the renormalization 
point. As far as we do not need to deal with the anomalous dimension of the A 2 operator, we can 
factorise the tree-level ghost propagator. The ghost dressing function is then written as: 

F R (q 2 ^ 2 ) = F R , peit (q 2 ,^) ( l + j /^Nl-l) ) ' (29) 

where the multiplicative correction to the purely perturbative Fn jPeit is determined up to correc- 
tions of the order 1/g 4 or hnq/fi (the Wilson coefficient at the leading logarithm is computed in 
appendix 03]) . 

5 The operators with an odd number of fields (d = 1,3/2; d^A and d^c) cannot satisfy colour and Lorentz 
invariance and do not contribute a non-zero non-perturbative expectation value, and cc does not contribute either 
because of the particular tensorial structure of the ghost-gluon vertex. 
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We can handle in the same way (see refs. [7, 8]) the OPE power correction to the gluon 
propagator and obtain 



w 




2 x 




3g 2 ,^(2) ,ab 



(G 



pert/ fiV 



(30) 



Then, after renormalization, one gets 



(G 
(G 



(2) 

_R,pcrt 



(A 2 ) 



Y b ( a 2 u 2 ) + (w ab ) v 

3 9 2 R (A 2 )Z 



+ 



(2) ) ab (a 2 a 



q 2 4(N 2 c - 1) 

and an appropriate projection gives for the gluon dressing function : 



(31) 



G R (q 2 ,n 2 ) = G R , pert (q 2 ,fi 2 ) 



3 9 2 R (A 2 ) K 
q 2 4(N 2 - 1 



(32) 



Finally, putting together the defining relation eq. (jHD and the results eqs. ( 12911321) we get 



«r(/^ 2 



lim ^F 2 (/i 2 ,A)G(/x 2 ,A) 



«r rt (9o) 



lim ^F 2 (q 2 ,A)G(q 2 ,A) F%(jjl 2 ,<&) G R (^ 2 ,q 2 ) 

A— >oo 47T 



a', 



(9o)-^l,pert(/ /2 ) So) ^fl.pert^ 2 , <?o) 1 



9 ffT(gQ 2 )Wi?, g g 
^ 4(iVg-l) 



(33) 



pert/ 2\ 



where q$ 3> Aqcd is some perturbative scale and the /3-function, and its coefficients in eq. ( 12TI) . 
of course describe the running of the perturbative part of the evolution, a^ n . 

The Wilson coefficient at the leading logarithm for the T-scheme MOM running coupling is 
presented in appendix [Bj where we also show that the inclusion of the logarithmic correction 
would induce no significant effect, provided that the coupling multiplying A 2 inside the bracket 
is taken to be renormalized also in T-scheme. Thus, for the sake of simplicity, eq. ( 1331) will be 
applied for our analysis in the next section. 
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3 Data Analysis 



In the following, we will first propose a "plateau" -procedure exploiting eq. fl33|) to get a reliable 
estimate of the Agc^-parameter from the lattice and we will apply it to previously published 
quenched lattice data [21,22] as a check of the method. 

3.1 The "plateau" method 

The goal being to get a trustworthy estimate of the A^g-parameter, one could attempt to do it by 
inverting the perturbative formula eq. (I23p and using in the inverted formula the lattice estimates 
of the running coupling obtained by means of eq. ([9]) for as many lattice momenta as possible. 
Then, one should look for a "plateau" of A^g in terms of momenta in the high-energy perturbative 
regime (this was done with the coupling defined by the three-gluon vertex in [4,5]). In the next 
subsection, fig. [21(a) shows the estimates of A^g so calculated for the lattice data presented in 
ref. [21,22] over 9 <p 2 < 33 GeV 2 . 

However, in order to take advantage of the largest possible momenta window one can use 
instead eq. ( 133]) . In this way we shall hopefully be able to extend towards low momenta the region 
over which to look for the best possible values of the gluon condensate and of A^jg 0. In other 
words, one requires the best-fit to a constant of 

fayi) = (p-,A(«i)) , 

with: a % = (34) 

1 + - 

Pl 

where A(a) is obtained by inverting the perturbative four-loop formula, eq. (12*31) . and c results 
from the best-fit (it appeared written in terms of the gluon condensate in eq. (I33|) ). Of course, 
A(a) reaches a "plateau" (if it does) behaving in terms of the momentum as a constant that we 
will take as our estimate of Aj^g. 

3.2 Applying the method 

The lattice data that we will exploit here to check the method we have explained above were 
previously presented in ref. [21]. We refer to this work for all the details concerning the lattice 
implementation: algorithms, action, Faddeev-Popov operator inversion, etc. 
The parameters of the whole set of simulations are described in table [T] 

6 This increases the statistics and reduces errors. It also avoids some possible systematic deviation appearing 
when lattice momentum components, in lattice units, approach 7r/2 (Brillouin's region border). 
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3 


Volume 


a -1 (GeV) 

\Ai \ vJ V> V / 


Number of confs 


6.0 


16 4 


1.96 


1000 


6.0 


24 4 


1.96 


500 


6.2 


24 4 


2.75 


500 


6.4 


32 4 


3.66 


250 



Table 1: Run parameters of the exploited data [21]. 



3.2.1 The scaling from different lattices 

It should first be noted that the scaling of eq. ([H]) from the several lattices we use is indeed 
satisfactory. The prescription of taking the infinite cut-off limit in eq. @ means in practice to 
have the lattice artifacts under control. This is in fact the case for UV ones. In particular, the 
hypercubic artifacts behaving as 0(a k YliVi) f° r the lattice propagators we analyze were cured, as 
explained in [21], by exploiting the i7 4 -symmetry 

As an indirect way of testing that scaling, we consider all the lattice propagators as functions of 
the momentum measured in lattice units, {i.e. with dimensionless momenta pLat — where 
a((3) is the lattice spacing in physical units at the particular bare lattice coupling g 2 = 6/6), 
and determine the ratios of a(/9)'s for the scaling to work. Then, still working in lattice units, 
the best-fit parameters to be obtained by applying the "plateau" -method will be a((3) A^r and 
a 2 (/3)g 2 ^(A 2 ) ^, and the ratio of those best-fit parameters for different lattices will provide the ratio 
of the corresponding lattice spacings. 



p 


Volume 


a(3)/a(6.2) (this work) 


a(P)/a(6.2) [1] 


deviations (%) 


6.0 


16 4 


1.368 


1.378 


0.7 


6.0 


24 4 


1.322 


1.378 


4.1 


6.2 


24 4 


1 


1 





6.4 


32 4 


0.768 


0.751 


2.2 



Table 2: Comparison of lattice spacings ratios obtained by means of the scaling of eq. © as explained 
in the text and of the string-tension method. 

In tab. [21 the ratio of lattice spacings obtained by the standard string-tension method [1] 
are compared with those obtained as explained above. More precisely : (i) we first determine 
Ajyfg-a(6.2) and a 2 (6.2) (A 2 ) r for the lattice data with (3 = 6.2; (ii) then, for each new (3, we 
determine x = a(6) / a(6.2) in such a way that a "plateau" for x a(6.2) Aj^ is obtained with a gluon 
condensate given by x 2 a 2 (6.2) g^(A 2 )ji. They agree very well, at least for the ratios computed for 
the three lattice simulations with roughly the same physical volume: (3 = 6.0(L = 16 = 1.58 fm), 
(3 = 6.2(L = 24 = 1.72 fm), (3 = 6.4(L = 32 = 1.72 fm). A slightly larger discrepancy (~ 4%) 
appears when comparing with data for the largest lattice ((3 = 6.0, L = 24 = 2.37 fm). We suspect 
that this is the manifestation of a finite- volume effect. Actually, if we compare the two simulations 
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at (3 = 6.0 for different volumes (see fig. []]), such an effect can be seen, although it decreases as 
the physical momentum increases (and becomes in practice negligible at p 2 ~ 9 GeV 2 ). 

Thus, one can conclude that the scaling of the coupling defined by eq. (0) for p 2 >9 GeV 2 is 
very good. Conversely, this argument provides an alternative method to determine the lattice size 
for a simulation at a given f3 in terms of the one known in physical units at any other one. 
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Figure 1: (a) Plot of ay defined by eq. (|9|) in terms of the square of the renormalization momentum as 
computed from the two lattices at /3 = 6.0 with different volumes: V = 16 4 (green boxes) and V = 24 4 
(red boxes), (b) A zoom onto the high momenta region of the left plot. 



3.2.2 Looking for the "plateau" 

In fig. El (a), we show the estimates of A^jg obtained when interpreting the lattice coupling com- 
puted by eq. © for any momentum 9 <p 2 < 33 GeV 2 in terms of the inverted four- loop perturbative 
formula for the coupling, eq. f[2"3l . The estimates systematically decrease as the squared momen- 
tum increases until around 22 GeV 2 ; above this value, only a noisy pattern results. In fig. EJ(b), 
the same is plotted but inverting instead the non-perturbative formula including power correc- 
tions, eq. ( 1331) . The value of the gluon condensate has been determined by requiring a "plateau" 
to exist (as explained in the previous section) over the total momenta window. 

One should realize that, had we not taken into account the noisy ballpark of points above 22 
GeV 2 and had we considered the perturbative regime as reached at that momentum, we would 
have got an estimate of A^jg roughly 35-40 MeV above the one obtained from the non-perturbative 
formula. In other words, the non-perturbative analysis seems to indicate that the perturbative 
regime is far from being achieved at p = 5 GeV. This is illustrated in figure [3] in which, adopting 
for Ajyjg the value 224 MeV which results from the non-perturbative analysis, we plot against the 
square of the renormalization momentum the coupling constant as computed by means of the non- 
perturbative formula ( 1331 (red curve) and of the perturbative one (123]) (blue curve). Displayed 
are also the lattice data, i.e. the values of ar obtained from eq. In figure [31 a the range in 
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(a) (b) 

Figure 2: (a) Plot of A^g (in GeV) computed by the inversion of the four-loop perturbative formula 
eq. (I23p as a function of the square of the momentum (in GeV 2 ); the coupling is estimated from the lattice 
data through eq. ([9]). (b) Same as plot (a) except for applying the non-perturbative formula eq. ([33]) for 
the coupling and looking for the gluon condensate generating the best plateau over 9 < p 2 < 33 GeV 2 . 

/i 2 one sees that the non-perturbative approach provides a fairly good agreement with the data, 
the x 2 being 1.3 per degree of freedom. On the contrary there is a clear disagreement with the 
perturbative formula. Furthermore, one can extrapolate the value of the up to very high 
momenta with eq. ( l33l . p 2 ~ 300 — 500 GeV 2 , where the purely perturbative eq. (123]) and the 
non-perturbative eq. ( 133]) . both with the same Ajjjg, generate in practice the same results. The plot 
of fig. [3j(b) shows indeed that the curve for the coupling extrapolated in this way joins perfectly 
the lattice estimates at high momenta taken from [9]. Thus, the inclusion of the non-perturbative 
OPE power correction, eq. (1331) . to describe the running of the coupling eliminates effectively the 
observed systematic deviations for the estimates of A^g from the momenta window from 3 GeV to 
5 GeV (fig. El (a) ) and essentially leads to the same estimate as was found from the perturbative 
regime at very high momentum. 

Thus, we have been able to obtain simultaneous best-fit values for both the gluon condensate 
and Ajyfg. It is however manifest that they are correlated by their determination: the larger the 
gluon condensate is, the smaller the value of A^jg has to be. In fig. HI we plot the ellipsoid defined 
by El x 2 (Ams, g^{A 2 )^) = xLm + 1 for a fitting window defined by p 2 > 8 GeV 2 and for one 
restricted to p 2 > 14 GeV 2 . It is seen that, neglecting other sources of errors like, for instance, 
the calibration of the lattices, but being conservative with the choice of the fitting window, one 

7 The errors on the lattice estimates of the coupling that were used to compute \ 2 were obtained by propagating 
the ones computed through the jackknife method for F and G in [21]. 
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Figure 3: (a) Plot of ay defined by eq. ([9]) in terms of the square of the renormalization momentum: 
the red solid line is computed with eq. ([33]) with A^jg = 224 MeV, the blue one with eq. (f23j) for the 
same A^jg and the data are obtained from the lattice data set-up in table [TJ (b) The same but with some 
additional lattice estimates for the coupling at very high momenta (300-500 GeV 2 ) taken from [9]. 

can conclude that our best-fit parameters incorporating only statistical errors are: 

Ag =0 = 224+^MeV 
g 2 T (A 2 ) R = 5.1+?;! GeV 2 . (35) 

These values are in very good agreement with the previous estimates from quenched lattice sim- 
ulations of the three-gluon Green function [7, 8] or, in the case of A^jg, from the implementation 
of the Schrodinger functional method [2] , although slightly larger than the one obtained by the 
ratio of ghost and gluon dressing functions [22] (see fig. Hl(b) and tab. [3]). Concerning the gluon 
condensate estimate only, it is worth pointing that it can be computed at the renormalization 
momentum /i 2 = 100 GeV 2 (see tab. [3]) and it also agrees very well with the estimate from the 
analysis of the quark propagator vector part, Z^, that gives: (A 2 )n :fl=w GeV = 1.76(8) GeV [26]. 

As a final remark, had we taken into account the leading-logarithm behaviour of the Wilson 
coefficient for the running coupling (applied eq. (PHI) ) instead of eq. (133]) ), the parameters so 
fitted would not significantly differ from those in eq. (1351) : we estimate a difference of ~ 4% in the 
determination of g^(A 2 )f( and less than 0.5% in that of A^g. 



4 Conclusions 

In the present paper we reconsider in some detail the determination of A^g from gluon and ghost 
Green functions using the MOM scheme. We stick here to the quenched case, or rather to the pure 

8 we define the errors by taking the larger ellipsoid and this could be maybe considered as to give account of 
some systematic effect related to the choice of the fitting window. 
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Average 




Figure 4: (a) The ellipsoid defined by x 2 (A^jg, g\ (A 2 )r) = x^in + 1- The y-axis is for expressed 
in GeV and x-axis for g^(A 2 )ji in GeV 2 . The small ellipsoid is obtained for a fitting window defined by 
p 2 > 9 GeV 2 and the larger is for p 2 > 14 GeV 2 . (b) Comparison with previous estimates of Aj^g in pure 
Yang-Mills collected in tab. [3j the blue triangle stands for the estimate in this work and the red square 
for the average of the others. The l-er error interval for the average (dashed red line) were estimated by 
treating the errors in tab. [3] as purely statistical ones. 

Yang-Mills SU (3) theory, having of course in mind to apply what we learn also to the unquenched 
situation. 

4.1 ghost-gluon vertex 

We give some details about the proper renormalisation of the ghost-gluon vertex in the MOM 
scheme mainly because we realised that there is some carelessness in literature. An obvious remark 
is that applying MOM to a vertex function needs to specify the kinematics of the renormalisation 
point. Renormalising at the scale [i may be performed in the symmetric case, with the three 
momenta at the renormalisation scale (p 2 = /i 2 ) or in the soft gluon limite (p g i UO n = ®,Pghost = A* 2 )' 





F 2 G (this work) 


Asym. 3-g [8] 


Sym. 3-g [8] 


F/G [22] 


[2] 


Ams (MeV) 


224^ 


260(18) 


233(28) 


270(30) 


238(19) 


^{A 2 ) R ^ (GeV) 


1.64(17) 


2.3(6) 


1.9(3) 


1.3(4) 





Table 3: Comparison of estimates of A^jg obtained from the analysis of the ghost-gluon vertex in this 
work (first column), the asymmetric 3-gluon vertex (second), the symmetric 3-gluon vertex (third), the 
ratio of gluon and ghost dressing functions (fourth) and with the Schrodinger functional method (last). 
The gluon condensate (A 2 )ji tfl has been obtained at the renormalization momentum /i = 10 GeV, for the 
sake of comparison with the other estimates, from eq. (|35p by applying g 2 (/d 2 = 100 GeV 2 )/4-7r = 0.15. 
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or with a vanishing incoming ghost momentum, etc. The latter case is the one in which Taylor's 
theorem applies which leads to Z\ — 1. We present in section I2TT1 an alternative derivation of the 
perturbative renormalisation of the coupling constant in the latter scheme, defined by eq. eq. (P), 
in agreement with the result by Chetyrkin [18]. The other kinematics lead to a finite but non trivial 
Z\ = 1 + 0(a 2 ). This difference has been often overlooked, presumably because it is assumed to 
be small. However, as we have shown in section 12.11 applying Z\ = 1 to the symmetric case leads 
to a 14 % systematic error on A^jg while it gives 23 % when applied to the soft gluon limit. 

4.2 The Ajjg plateau 

Ajyfg is a constant independent on the scale /i. Inverting the perturbative expansion of the coupling 
constant one can invert eq. eq. (l23l) leading for each /i to A^jg(/i 2 ) from a<r(/i 2 ) B If we were in a 
perturbative region of \x Aj^g(/x 2 ) should not depend on \i up to statistical errors. One should see 
a nice "plateau". Fig. El (a) shows that this is far from being the case up to /i 2 = 30Ge^ 2 . We 
have since long advocated that there is a sizeable non-perturbative contribution from the vev of 
the unique (in Landau gauge) dimension 2 operator (A 2 ). We propose to fit this condensate by 
adjusting the resulting A^jg to a " plateau" . This is successfully achieved, see fig. El(b). Since we 
scan a large window in the scale \i we believe that we are in a position to claim that we indeed see 
a non-perturbative 0(1/ n 2 ) contribution rather than the effect of logarithmically behaved higher 
orders in perturbation theory (0(a 5 )). 

4.3 Comparison of different estimates of Ajj$ 

We have performed a comparison of different estimates of A^ and (A 2 ) in the pure Yang-Mills 
theory using the coupling constant defined in eq. (151) , the MOM coupling constant from symmetric 
three gluon vertex function, the MOM coupling constant from the three gluon vertex function with 
one vanishing momentum and from the ghost to gluon propagator ratio, and also with the estimate 
of Aj^g from the Schrodinger functional approach. The result is reported in table [3] and fig. Hl(b). 
The agreement is quite satisfactory. Fig [31(b) shows also a good agreement of our fit from ay^ 2 ) 
with very large \x measurements from [9]. Notice also that A^g from a^(/i 2 ) has the smallest 
statistical errors due to the fact that it relies only on propagator, not on noisier three point Green 
functions. 

This opens a possibility of using the matching of Aj^g as computed from different lattices 
in order to fit the lattice spacing ratio. One might also match directly ar(/i 2 ) from different 
lattices, a procedure which is not constrained to large scales and does not need to estimate the 
(A 2 ) condensate. In fact from eq. Q we get directly a quantity which should be independant 
of the lattice spacing at the same \x in physical units, up to 0(l/a 2 ) artifacts. This method is 
complementary to the use of Sommer's parameter ro [27] and it also only depends on gauge fields. 

9 This can be done in any MOM scheme using the apropriate equivalent to eq. (|23jl . 
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A Appendix: ghost and gluon propagators anomalous di- 
mension in MOM 



The ghost and gluon anomalous dimension can be computed in MOM scheme by applying eq. (IT7|) 
with the results obtained in MS for the radiative corrections of all the relevant Green functions [18, 
20]. Thus, one obtains for the coefficients defined in eq. ( 1T61) : 



To 
7i 



73 



9 
4 

813 _ l3N f 
Iff 4~ 



157303 14909iV7 125N] 5697C(3) 21 
72 _ 64 48 + "Ts 32 T f ^ ' 



219384137 30925009iV f 288155iV| 2705iVj 9207729C(3) 
1536 1152 + 216 162 512 

+ ^ /C( 3) - ^ m - + 1^5 W(B) (36) 



70 = Y--3 



7i 



73 



13 2N f 

2 



3727 250 Nf 20 Nj 
~24 9 + 27 
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288 16 V 3 3 / 

2 /1681 16 C(3) \ 200 N} 
7 V 18 9 J 243 

3011547563 _ 18987543C(3) 1431945C(5) 

6912 256 64 

^ N / 221198219 2897113C(3) 845275C(5)^ 



-iVj 



1728 + 216 + 96 J 
'6816713 60427C(3) 4640C(5)\ 



648 

373823 




) 



+ Ar H-w^ l+ — (37) 

These coefficients appear for the expansion, given by eq. (1171) . of the MOM-renormalized ghost and 
gluon anomalous dimension in terms of the MS-coupling. However, provided that the /3-function 
for any other renormalization scheme is known, it can be applied to replace a^g in eq. fflBI) by 
the coupling in that scheme. 
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B Appendix: Wilson coefficients at leading logarithms 

The purpose of this appendix is to present up to leading logarithms the subleading Wilson coeffi- 
cients in eqs. f l2TIll3"2"l) and, in view of checking the validity of neglecting those logarithms, estimate 
their impact on the momenta window we use for our fits. Following [8], let us write 

F R (q 2 , /i 2 ) = co cV) j + c 2 (j^, J (* 

for gluon and ghost propagators. Then, with the help of the appropriate renormalization constants 
one can rewrite eq. (1381) in terms of bare quantities: 



q 2 


a(/i 2 )^ 




V 2 


A(N* - l)q 2 


q 2 


a(/i 2 )^ 




V 2 


A(N 2 C - l)q 2 



(q 2 ,A 2 ) = Z 3 (fi 2 ,A 2 )c [^,a(^ 



+ Z^ 2 ,A 2 )Z-J(^,A 2 )c 2 (^^ 2 ?) m fl\ )q2 , (39) 



where A\ = ZJ.A 2 . A totally analogous equation for the ghost dressing function F(q 2 ,A 2 ), with 
q and Z 3 in place of q and Z 3 . Now, as the /^-dependence of both l.h.s. and r.h.s. of eq. f[3"9"]) 
should match each other for any q, one can take the logarithmic derivative with respect to /1 and 
infinite cut-off limit, term by term, on r.h.s. and obtains: 

7(«(/i 2 )) + liTT^ + /9(«(/i 2 ))|-l \nc (^a^ 2 )) = 



k d log /1 2 da J \ /1 2 ' 

-7A 2 («(/i 2 ))+ 7 («(/i 2 )) + {^^ + /3(a(/x 2 ))^} lnc 2 ^,a(/i 2 )^) = 0, (40) 

where 7(a(yU. 2 )) is the gluon propagator anomalous dimension defined in eq. ( 1T61) and 

7^(<V)) = lim -^ln^ 2 (/i 2 ,A 2 ) = - 7( f 2 ^fl + ... (41) 
a->oo d In fi z Air 

Both eqs. (1401) can be finally combined to give: 



- 7 ^(<V)) + — — — ^ - + (3(a(fi 2 ))-^— — ^ f = . (42) 



<91og/i 2 <9aJ Co ^ ;a(/i 2 



We can proceed in the same way for the ghost dressing function and derive analogous equations 
for the Wilson coefficients, Cj, that differ from those for q only because 7(a(/i 2 )) takes the place of 
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7(a(/U 2 )). Thus, the combination c^/cq obeys exactly the same eq. ( l42l) . above derived for C2/C0, 
that can be solved at the leading logarithm as explained in [8] to give: 



c 2 



a(fj, 2 



c 2 



("0 



a{\i 2 



Co 



a{\i 2 



g\g 2 



-7<f /A) 



(43) 



The boundary condition comes from requiring eq. (1401) to be equal to eq. ( 1291) for the ghost and 



eq. ( 1321) for the gluon at /x 2 = g 2 . The coefficient 7^ was computed to be 35/4 for the first 
time in [8]. Of course, eqs. (T4T)1) define not only the dependence of the Wilson coefficient on the 
renormalization momentum, /i 2 , but also that on the momentum scale q 2 because of standard di- 
mensional arguments: the only dimensionless quantities^) are the ratio q 2 / fi 2 and a. Then, putting 
all toghether, the non-perturbative formula for the running coupling at the leading logarithm is 
given by 



a 



pert 



(M 2 



/ 

1 - 




QCD 


-9/44 

1 


V 




\ QCD / 





g 2 M)(A 2 ) RA 
4(iV 2 - 1 ^ 



\ 



(44) 



/ 



where the only correction to eq. (!33l comes from the ratio of logarithms inside the bracket that, 
as can be seen in fig. El introduces no significant deviation. 




Figure 5: a NP /a pert in terms of the square of the momentum computed by using both eq. 
blue) and eq. ([33]) (solid red). 



(dashed 



10 Other dimensionless quantities can be obtained with the help of Aqcd, but this is a non-perturbative parameter 
not emerging in the Wilson coefficient dominated by the short-distance singularities of the OPE expansion and 
only coding pcrturbative information in the SVZ approach. 
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